Frustrated square lattice with spatial anisotropy: 
crystal structure and magnetic properties of PbZnVO(P04)2 
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Crystal structure and magnetic properties of the layered vanadium phosphate PbZnVO(PC>4)2 
are studied using x-ray powder diffraction, magnetization and specific heat measurements, as well as 
band structure calculations. The compound resembles AA'VO(P04)2 vanadium phosphates and fits 
to the extended frustrated square lattice model with the couplings Ji , J[ between nearest-neighbors 
and J2,J% between next-nearest-neighbors. The temperature dependence of the magnetization 
yields estimates of averaged nearest-neighbor and next-nearest-neighbor couplings, Ji ~ —5.2 K 
and J2 — 10.0 K, respectively. The effective frustration ratio a = J2/J1 amounts to —1.9 and 
suggests columnar antiferromagnetic ordering in PbZnVO(P04)2. Specific heat data support the 
estimates of Ji and J2 and indicate a likely magnetic ordering transition at 3.9 K. However, the 
averaged couplings underestimate the saturation field, thus pointing to the spatial anisotropy of 
the nearest-neighbor interactions. Band structure calculations confirm the identification of ferro- 
magnetic Ji,J{ and antiferromagnetic 3i,J'i in PbZnVO(P04)2 and yield (J[ — Ji) ~ 1.1 K in 
excellent agreement with the experimental value of 1.1 K, deduced from the difference between the 
expected and experimentally measured saturation fields. Based on the comparison of layered vana- 
dium phosphates with different metal cations, we show that a moderate spatial anisotropy of the 
frustrated square lattice has minor influence on the thermodynamic properties of the model. We 
discuss relevant geometrical parameters, controlling the exchange interactions in these compounds, 
and propose a new route towards strongly frustrated square lattice materials. 

PACS numbers: 75.50.-y, 75.40.Cx, 75.30.Et, 75.10.Jm 



I. INTRODUCTION 

Quantum magnetism is one of the active research 
topics in solid state physics. Quantum effects lead to 
numerous unusual properties, interesting with respect 
to the fundamental understanding of critical phenom- 
ena (spin-liquid ground states, Bose-Einstcin condensa- 
tion in high magnetic fields)^ and potential techno- 
logical applications (magnetoelectric coupling, ballistic 
heat transport) The search for new quantum magnets 
that enable to challenge theoretical predictions remains a 
long-standing problem in solid state science, because the- 
oretical studies usually deal with relatively simple mod- 
els, like the Heisenberg model, while the properties of 
real compounds are often determined by the interplay 
of numerous effects, such as isotropic and anisotropic 
exchange interactions, orbital and/or lattice degrees of 
freedom, etc. 

The frustrated square lattice (FSL) with isotropic ex- 
change is an archetypal spin model in two dimensions. 5 
This model reveals strong quantum fluctuations due to 
the combination of low-dimensionality and frustration. 
The specific regime of the model is determined by the 
frustration ratio, i.e., by the ratio of the next-nearest- 
neighbor (NNN) interaction (J 2 ) to the nearest-neighbor 
(NN) interaction (Ji): a = J2/ J\- The frustration ra- 
tio determines the ground state of the system - either- 
ordered [Neel antiferromagnetic, columnar antiferromag- 
netic, or ferromagnetic (FM) ordering] or disordered.®' 7 



The precise nature of the disordered [presumably spin- 
liquid at J2/J1 ~ 0.5 (Ref. H) and nematic at J2/J1 — 
—0.5 (Ref. @)] ground state remains controversial, be- 
cause a theoretical treatment of the model is inevitably 
approximate. 

Experimental studies aim at the search for materials 
that fit to the FSL model. Recently, extensive inves- 
tigations identified a number of vanadium-based FSL 
compounds. Li 2 VOSi04 and Li 2 VOGeC>4 reveal anti- 
ferromagnetic (AFM) Ji and J 2 with J 2 » J X ^ML 
Thus, the columnar AFM ordering is established, while 
the frustration is weaki^ 2 . In the J\ > J 2 (both J\ 
and J 2 AFM) region, VOM0O4 undergoes Neel-type 
ordering ; 13 i 14 yet the presence of the long-range order 
in PbV0 3 is still under debate,! 5 -^ The region of FM J Y 
~ AFM J 2 is probed by the family of layered vanadium 
phosphates AA'VO(P0 4 ) 2 (AA' = Pb 2 , BaZn, SrZn, and 
BaCd) i 17 ' 18 i 19 ' 20 None of the latter compounds fall in the 
critical region with J%/ J\ ~ —0.5, although two of them, 
BaCdVO(P0 4 ) 2 and SrZnVO(P0 4 ) 2 , reveal a ~ -1 and 
show a pronounced effect of the frustration on the ther- 
modynamic properties ! 18 ! 19 An additional feature of the 
layered vanadium phosphates is the lack of the tetrago- 
nal symmetry and the resulting spatial anisotropy of the 
spin lattice (both for the NN and NNN couplings) i2I 

Motivated by the above-mentioned demand for 
strongly frustrated FSL materials, we attempted to ex- 
tend the family of layered vanadium phosphates. Metal 
cations accommodated between the V-P-0 layers do not 
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take part in the magnetic exchange. However, the size of 
these cations determines the geometry of the magnetic 
layer,— hence the replacement of metal cations can be 
a promising way to tune the spin system of the mate- 
rial. The approach of cation substitution is widely used 
in solid state chemistry. For example, the replacement 
of Li by Na in L1VO2 leads to a change in the orbital 
ordering pattern and results in the long-range magnetic 
ordering in NaVC>2 (Ref. [22] ) instead of the formation 
of trimer clusters in L1VO2. 23 Following a similar ap- 
proach for the FSL compounds, we succeeded in prepar- 
ing the vanadium phosphate PbZnVO(P04)2 that be- 
longs to the family of layered AA'VO(P04)2 vanadium 
phosphates and combines certain structural features of 
SrZnVO(P0 4 ) 2 and BaZnVO(P0 4 ) 2 . To the best of 
our knowledge, neither the preparation and the crystal 
structure, nor the magnetic properties of PbZnVO(P04)2 
have been reported before. The only exception is the 
high-field magnetization curve of PbZnVO(P0 4 )2 pre- 
sented in the comparative study of high-field properties 
of the FSL compounds.— 

In the following, we present the results of a combined 
- structural, phenomenological, and microscopic - study 
of the model FSL compound PbZnVO(P0 4 ) 2 . We start 
with the methodological aspects in Sec.|TTJ In Sec. lIIIl the 
structural data are reported. Sec. IIVI deals with the ex- 
perimental study of the magnetic properties, while Sec.lVl 
presents band structure calculations and the evaluation 
of individual exchange couplings. In Sec. IVI1 we fur- 
ther exploit the computational approach and consider the 
influence of different structural factors on the magnetic 
interactions in layered vanadium phosphates. Sec. IVIII 
concludes the study with the discussion and a summary. 



II. METHODS 

Polycrystalline samples of PbZnVO(P04)2 were ob- 
tained by heating a mixture of PbZnP207, V2O3, and 
V2O5 in an evacuated quartz tube (10~ 2 mbar) at 700 °C 
for 24 hours. Phase composition of the samples were 
controlled using x-ray diffraction (Huber G670 Guinier 
camera, CuK Q i radiation, 29 = 3 — 100° angle range). 
PbZnP207 was obtained by heating a stoichiometric mix- 
ture of PbO, ZnO, and NH 4 H 2 P0 4 in air at 750 °C for 
48 hours. 

The best sample was obtained from the reactant 
mixture with the stoichiometric cation composition 
and the slight oxygen excess corresponding to the 
PbZnVOi.o9(P0 4 )2 formula. The sample contained the 
targeted PbZnVO(P0 4 )2 phase and the minor impurity 
of diamagnetic PbZnP207 (about 2 wt. %). This sam- 
ple was further used for the structure refinement and 
for the thermodynamic measurements. According to the 
structure refinement and magnetization measurements, 
the PbZnVO(P0 4 )2 compound is stoichiometric^ The 
annealing of the PbZnVO(P0 4 )2 samples above 700 °C 
resulted in the decomposition of the compound. 



The x-ray powder pattern for the structure refine- 
ment was collected using the STOE STADI-P diffrac- 
tometer (transmission geometry, Ge(lll) monochroma- 
tor, CuK Q i radiation, linear position-sensitive detector, 
angle range 29 = 7 — 100°). The structure refinement 
was performed using the JANA2000 program. 26 

Temperature dependence of the magnetization was 
measured using the Quantum Design MPMS SQUID in 
the temperature range of 2 — 360 K. Heat capacity data 
were collected with the commercial PPMS setup. To 
improve the contacts between the grains, the pellet for 
the heat capacity measurement was additionally annealed 
overnight in the evacuated quartz tube at 500 °C. The 
x-ray study of the re-annealed sample evidenced that the 
phase composition remained unchanged. 

Scalar-relativistic band structure calculations were 
performed using the full-potential local-orbital scheme 
(FPL08. 50-32) within the local density approximation 
(LDA) of density functional theory^ The Perdew- 
Wang version of the exchange-correlation potential was 
applied-S. The calculations were performed for the crys- 
tallographic unit cell of PbZnVO(P0 4 ) 2 (112 atoms, k 
mesh of 256 points with 75 points in the irreducible 
part) and for a number of model structures with the 
formal composition of LiVOP0 4 (32 atoms, k mesh of 
4096 points with 729 points in the irreducible part). The 
construction of the model structures is further described 
in Ref. [2l|. LDA band structures were used to select 
the relevant states and to analyze the respective bands 
with tight-binding (TB) models. The parameters of these 
models were evaluated as overlap integrals of Wannicr 
functions localized on vanadium site a 29 ' 30 and further 
used to estimate exchange integrals. Details of the pro- 
cedure are given in Sec. |Vj 



III. CRYSTAL STRUCTURE 

The X-ray powder pattern of PbZnVO(P0 4 )2 was 
indexed in an orthorhombic unit cell with lattice 
parameters a = 8.763(1) A, b = 9.072(1) A, 
c = 18.070(5) A. This unit cell is similar to unit cells of 
the compositionally-related compounds BaZnVO(P0 4 )2 
(a = 8.814 A, b = 9.039 A, c = 18.538 A)2I 
and SrZnVO(P0 4 ) 2 (a = 9.066 A, b = 9.012 A, 
c = 17.513 A). 32 Reflection conditions (hkO), h = 2n; 
(h0l),l = 2n;(0kl),k = 2n unambiguously pointed to 
the Pbca space group and supported the similarity to 
the AA'VO(P0 4 )2 layered vanadium phosphates (AA' 
= BaZn, SrZn, BaCd). To select the proper starting 
model for the structure refinement, we simulated pow- 
der patterns using lattice parameters of PbZnVO(P0 4 )2, 
atomic positions of BaZnVO(P0 4 ) 2 or SrZnVO(P0 4 ) 2 , 
and Pb atom instead of Ba/Sr. The structural data 
for BaZnVO(P0 4 )2 showed better agreement with the 
experimental pattern. Therefore, the structure of 
BaZnVO(P0 4 )2 was used as a starting model for the re- 
finement. 
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TABLE I: Atomic positions and isotropic atomic displacement 
parameters (U iso , in units of 10" 2 A 2 ) for PbZnVO(P0 4 ) 2 . 
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FIG. 1: (Color online) Structure refinement for 
PbZnVO(P04)2: experimental, calculated, and differ- 
ence X-ray diffraction patterns. Upper and lower sets of 
ticks show line positions for the PbZnP207 impurity and for 
PbZnVO(P0 4 ) 2 , respectively. 



For the final refinement, the PDZ11P2O7 impurity was 
taken into account: 33 Atomic displacement parameters 
for two positions of phosphorous and for nine positions 
of oxygen were independently constrained due to the 
huge difference in the scattering powers of Pb and the 
light O and P atoms. The refinement converged to 
Ri = 0.019, R p = 0.028, and X 2 = 2.55. The experi- 
mental, simulated, and difference diffraction patterns are 
shown in Fig. [T] Atomic positions are listed in Table Q] 
while interatomic distances, relevant for a further discus- 
sion of the magnetic interactions, are presented in Ta- 
ble nn 

The crystal structure of PbZnVO(PC>4)2 is shown in 
Fig. [21 Vanadium atoms form distorted octahedra with 
the short vanadyl bond of 1.55(1) A, four longer bonds 
of 2.00 — 2.13 A in the equatorial plane, and a rather 
long bond of 2.38(1) A opposite to the vanadyl one. If 
the latter atom is cast away, the vanadium coordina- 
tion is reduced to a square pyramid (in the following, we 
will refer to the square-pyramidal coordination of vana- 
dium in order to simplify the comparison to the other 
AA'VO(P04)2 compounds). Phosphorous atoms reside 
in slightly distorted PO4 tetrahedra typical for phosphate 
compounds. The lead cation has seven oxygen neighbors 
with Pb-0 distances ranging from 2.47 A to 2.99 A. The 
wide range of the Pb-0 distances indicates the asymme- 
try of the local environment due to the presence of the 
6s 2 lone pair. Finally, zinc atoms form nearly regular 
tetrahedra with an additional oxygen atom at a longer 
distance of 2.38(1) A. 

Similar to the other AA'VO(P04)2 compounds, there 
are two incquivalcnt phosphorous sites in the crystal 
structure. The P(1)C>4 tetrahedra share corners to 
the VO5 pyramids and form [VOPO4] layers. The 
P(2)C>4 tetrahedra along with the Pb and Zn cations 
reside in the interlayer space and can be considered 
as two-dimensional [AA'P0 4 ] blocks (Fig. [2J|. The 



All the atoms occupy the general 8c position. 
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TABLE II: Relevant interatomic distances (in A) for the in- 
teractions J2 and J 2 . Vanadium is surrounded by six oxygen 
atoms. However, two of these atoms occupy axial apices of 
the octahedron and do not take part in the superexchange. 
Phosphorous is coordinated by four oxygen atoms, while the 
two O-O distances measure the edges of the PO4 tetrahedron. 
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[VOPO4] layers are buckled with the buckling angle 
Lp = 163° [compare to 160° in BaZnVO(P0 4 ) 2 and 150° 
in SrZnVO(PC>4)2]. The VO5 pyramids are connected 
via PO4 tetrahedra that give rise to four incquivalcnt 
exchange couplings: Ji,J( between nearest-neighbors 
and JiiJ'i between next- nearest- neighbors (left panel of 
Fig. f2J). Specific structural features of the magnetic layer 
along with the resulting exchange interactions constitute 
the main subject of the present paper and will be thor- 
oughly discussed in the following sections. In the re- 
mainder of this section, we will briefly comment on the 
structure of the [AA'PC^] blocks and on the role of the 
lead cation in PbZnVO(P04)2 as compared to the other 
AA'VO(P0 4 ) 2 compounds. 

The crystal structures of the AA'VO(PC>4)2 com- 
pounds contain A and A' cations that reside in cavities 
formed by the PO4 tetrahedra (right panel of Fig. EJ) • The 
arrangement of these tetrahedra is rather flexible and 
can be tuned for a specific metal cation. In the previous 
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FIG. 2: (Color online) Crystal structure and spin model of PbZnVO(P04)2: the [VOPO4] layer and the magnetic interactions 
(left panel); overall view of the crystal structure (middle panel); bonding of Pb and Zn cations in the interlayer blocks (right 
panel). Larger and smaller spheres denote Pb and Zn. The tp angle measures the buckling of the [VOPO4] layers. In the sketch 
of the spin lattice (left panel), solid, dash-dotted, dashed, and dotted lines show the interactions Ji, J[, J2, and J2, respectively. 



study, we have shown that barium coordinates all the four 
tctrahedra around the cavity, while strontium prefers a 
lower coordination number and links to three tetrahedra 
only (see Fig. 6 in Ref. HH). Within a naive picture, the 
lead-containing structure should resemble the strontium- 
containing counterpart due to the similar ionic radii of Pb 
(1.23 A) and Sr (1.26 k)M However, links to all the four 
tetrahedra around the cavity are present (right panel of 
Fig. [5]) . This feature explains the close similarity of the 
PbZnVO(P0 4 ) 2 and BaZnVO(P0 4 ) 2 structures, despite 
the smaller size of the Pb cation compared to Ba. The 
different behavior of lead and strontium is a well-known 
feature caused by the asymmetric local environment of 
lead (compare, e.g., SrVC>3 and PbVOs).— Yet the sim- 
ilarity between lead and barium is not common. 

We should also note that PbZnVO(PC>4) 2 belongs to 
the family of the AA'VO(P0 4 ) 2 phosphates with A ^ A'. 
The lead-zinc compound is rather different from the pure 
lead vanadylphosphate Pb 2 VO(P04) 2 with a monoclinic 
structured This indicates that the accommodation of 
the two asymmetric lead cations requires an overall dis- 
tortion of the crystal structure, while the single lead 
cation in PbZnVO(P04) 2 can be tolerated by simple 
tilts of the PO4 tetrahedra in the [AA'PC^] interlayer 
blocks. Indeed, the PO4 tetrahedra in SrZnVO(PC>4) 2 
and BaCdVO(P0 4 ) 2 have faces parallel to the [VOP0 4 ] 
layers, while this is not the case for the PbZnVO(PC>4) 2 
structure (compare the middle panel of Fig. [2] with Fig. 2 
in Ref. [HI). 



IV. THERMODYNAMIC PROPERTIES 

A. Magnetic susceptibility 

Magnetic susceptibility (x) vs. temperature (T) data 
for PbZnVO(PC>4) 2 reveal a maximum at about 8.8 K 
(see Fig. [3]). Such maxima are typical for low-dimensional 
spin systems and indicate the onset of short-range spin 



correlations. Below the maximum, field-dependent be- 
havior is observed. The data measured in low fields (be- 
low 1 T) show an anomaly at 3.9 K, while the data col- 
lected at 5 T present a bend at 4.2 K. These features 
are related to a phase transition that is also evidenced 
by a peak in the specific heat data (Fig. [5]). The field 
dependence of the transition temperature points to the 
magnetic origin of this transition. Based on the similarity 
to the other AA'VO(P04) 2 compounds, we can suggest 
that PbZnVO(PC>4) 2 undergoes antiferromagnetic order- 
ing at TV = 3.9 K. The transition temperature is slightly 
increased in the magnetic field due to the suppression of 
quantum fluctuations. A similar effect has been recently 
observed in BaCdVO(P0 4 ) 2 (Ref. (H) and in the un- 
frustrated square-lattice compound [Cu(HF 2 )(pyz) 2 ]BF 4 
(Ref. [37]). In the latter case, theoretical treatment was 
also given. Below the transition, the susceptibility re- 
mains nearly temperature-independent, as expected for 
a regular antiferromagnet. 

The high-temperature part of the susceptibility curves 
shows a Curie- Weiss (CW) behavior. Above 40 K, we fit 
the data with the expression 



c 



T + 6 C 



w 



(1) 



where the temperature-independent term Xo W accounts 
for the diamagnetic and Van Vleck contributions, C is 
the Curie constant, and 6*cw is the Curie- Weiss tem- 
perature. The fit (see the inset in the left panel of 
Fig. [3]) resulted in X q W = —1-8(1) • 1CT 4 emu/mol, 
C* = 0.354(1) emu K/mol, and 8 GW = 6.8(1) K. The 
C value corresponds to the effective magnetic moment 
of 1.64(1) /ib which is slightly reduced as compared to 
the expected spin-only value of 1.73 [is for V +4 . This 
reduction can be attributed to a slight reduction of the 
g-factor of vanadium due to the weak spin-orbit coupling, 
to the effect of the non-magnetic impurity phase, and to 
a very small amount of V +5 defects in PbZnVO(P04) 2 . 
To evaluate exchange couplings in PbZnVO(P04) 2 , 
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FIG. 3: (Color online) Temperature dependence of the magnetic susceptibility of PbZnVO(P04)2. Left panel: curves measured 
at 0.5 T (open circles), 1 T (open squares), and 5 T (filled circles); the inset shows the Curie- Weiss (CW) fit above 40 K. Right 
panel: HTSE fits (a') and (a") (solid lines), (&') and (6") (dashed lines) according to Table Hill the inset shows the difference 
between (a') and (a") in the high-temperature region, while dots indicate the experimental data. 



the susceptibility data are fitted with a high-temperature 
series expansion (HTSE) for the regular FSL model. 11 
In Sec. IIII1 we have shown that the structure of 
PbZnVO(P04)2 is rather complex, hence an extended 
spin model with four inequivalent couplings should be 
used. However, the spatial anisotropy (the difference be- 
tween J \ and J[ or J 2 and J 2 ) has a minor effect on ther- 
modynamic properties at sufficiently high temperatures, 
thus the expression for the regular model can be used. 
The resulting parameters should be considered as aver- 
aged couplings: J\ = (J 1 + J[)/2 and J 2 = {J 2 + J 2 )/2. 21 

The HTSE fits are known to be ambiguous. First, 
multiple solutions can be obtained, because individual 
couplings are separated in high orders of the expansion 
only, while the lowest, second-order term leads to a linear 
combination of the couplings (similar to the Curie- Weiss 
temperature #cw)- Second, the precise J values depend 
on the temperature range used. The temperature range 
of the fitting should lie within the convergence domain of 
the series, while the latter depends on the actual J values 
and makes the fitting procedure iterative. The ambigu- 
ity due to the multiple solutions can be handled via a 
reference to independent experimental data (saturation 
fielcU 9 - or specific heati£). Nevertheless, to remove the 
second ambiguity and to obtain accurate J values, one 
has to vary the lower boundary of the fitting range (T m ; n ) 
and to check the convergence of the series at T m j n . 

The experimental susceptibility data are readily fitted 
with the HTSE expression supplemented with an addi- 
tional temperature-independent xo term. Similar to the 
other AA'VO(P04)2 compounds ; 18 ' 19 we find two differ- 
ent solutions: (a) J\ < 0, J 2 > 0; and (b) J\ > 0, J 2 < 0. 
The J values for both solutions at different T m ; n are listed 
in Table ILLTI To discuss the dependence on T m ; n , we will 
use the solution (a) as an example. 

The value of T m ; n is varied between 8 K and 16 K. The 



upper panel of Fig. 0] shows that the J values are rather 
stable for T m ; n = 11 — 16 K. However, lower T m - m leads 
to a slight decrease of J 2 and to a notable increase of J\ , 
thus reducing the effective frustration ratio a = J 2 / 1 J\. 
To check the convergence of the series at T m i n , we use the 
ratio of the highest-order HTSE term to the susceptibility 
value. The change in the convergence is non-monotonous 
(see the middle panel of Fig. @}. As T min is decreased 
from 16 K to 11 K, the convergence is deteriorated. Yet 
the change in the frustration ratio for the solutions at 
low T min improves the convergence for T min = 9 — 10 K. 
Basically, we find two equally reliable solutions that cor- 
respond to T min = 16 K (a') and T min = 9 - 10 K (a"). 

Now, we will select the correct solution by calculating 
the saturation field and comparing it to the experimental 
value oi^H s = 23.4 At h/\J\\ > 0.5 [solution (a)], 
the saturation field of the regular FSL model is expressed 
as follows: 38 

= 2fcs(Ji + 2J 2 )/( gf i B ), (2) 

where the superscript "reg" denotes the saturation field 
for the regular model, and g is the ^-factor (we use 
g = 1.95 as a representative value for V +4 -containing 
compounds 3 - 9 ). Solution (b) corresponds to the Neel 
ordering, and the saturation field depends on J\ only: 
fioH r s cs — 4kB J\/{gnB)- The resulting [i H r s eg values are 
plotted in the bottom panel of Fig. 2J The solutions of 
type (6) overestimate the saturation field by 1.5 — 2.0 T 
and can be rejected, as further supported by band struc- 
ture results (see Sec. [V]) suggesting the FM Ji and AFM 
J 2 . The solution (a') underestimates the saturation field, 
while the solution (a") shows the best agreement with 
the experimental fioH s of 23.4 T. At first glance, this 
result suggests (a") as the final answer. However, one 
should be aware of the difference between the regular 
model assumed in ([2J and the real spin model with in- 
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TABLE III: HTSE fits of the magnetic susceptibility data. T m i n is the minimum temperature of the fitting range, xo is the 
temperature-independent contribution, g is the g-factor, Ji and J2 are averaged NN and NNN couplings, respectively, and 
fioHl ^ is the saturation field for the regular FSL model, as calculated via Eq. ((2j). (a) and (6) refer to the solutions with 
different signs of J\ and J2. The primes discriminate the solutions at different r m ; n . 



Fit 


T min (K) 


Xo (10 4 emu/mol) 


9 


Ji (K) 


7h (K) 


MoH s reg (T) 


(a') 


16 


-1.71(3) 


1.865(2) 


-5.2(2) 


10.0(1) 


22.6 


(a") 


10 


-2.11(4) 


1.889(2) 


-3.5(1) 


9.4(1) 


23.4 


(b') 


16 


-1.84(3) 


1.877(2) 


8.4(1) 


-2.6(1) 


25.6 


(6") 


10 


-2.08(4) 


1.892(2) 


8.3(1) 


-1.7(1) 


25.2 



equivalent couplings J% and J[ (Fig. [2]). In our previous 
workj2i we have shown that the difference between J\ 
and J[ increases I2qH s compared to noHl cg . Thus, both 
the solutions (a') and (a") are possible and would cor- 
respond to different regimes of the extended FSL model 
in PbZnVO(PC>4)2: (a') fits better to the square lattice 
with spatial anisotropy of the NN couplings, while (a") 
corresponds to the regular square lattice with J\ = J[. 

The final choice between (a') and (a") can be made af- 
ter a close examination of both fits. In the case of fit (a"), 
the low-temperature part of the data is fitted at the cost 
of a less accurate high-temperature behavior. Plotting 
the 1/x curves in the high-temperature region, we find 
that the (a") fit slightly deviates from the experimental 
data, while the fit (a') perfectly matches the data (see the 
inset in the right panel of Fig. |4j. This can also be seen 
from the xo values in Table IIIII compared to the Curie- 
Weiss fit value Xo W - -1-8 ■ 10" 4 emu/mol. The fit (a") 
overestimates xo, hence the high-temperature behavior of 
the experimental data cannot be accurately reproduced. 

The above arguments provide strong evidence for 
the solution (a 1 ) and for the scenario of the spatially 
anisotropic FSL. We further checked this conclusion by 
fitting the data collected in different fields. The J val- 
ues from different data sets coincided within 0.1 K, i.e., 
within the error bar of the fit. Thus, the analysis of 
the susceptibility data provides accurate estimates of J\ 
and Ji , while the further reference to the saturation field 
points to the spatial anisotropy of the NN couplings. 
According to Ref. [24], the difference between fioH s and 
fioHl cg is a measure of the spatial anisotropy: we find 
that J[ — J\ ~ 1.1 K. This value will be further com- 
pared with computational results in Sec. |Vj The satura- 
tion field could also be increased due to the AFM inter- 
layer couplings^ - However, our computational estimates 
suggest that the interlayer couplings are below 0.1 K (see 
Sec. |V| , hence their effect on the saturation field value 
will not exceed the experimental error bar of H s . 

Our analysis demonstrates that HTSE fits of the sus- 
ceptibility data yield accurate estimates of averaged ex- 
change couplings on the square lattice. Due to the sym- 
metry of the model hamiltonian, the fitting procedure 
is ambiguous. However, the reference to the experimen- 
tally measured saturation field can remove this ambigu- 
ity. Furthermore, it is crucial to collect the susceptibility 
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Minimal temperature of the fit: T mjn (K) 

FIG. 4: (Color online) Results of the susceptibility fit depend- 
ing on the minimum temperature of the fitting range (T m i n ). 
Upper panel: averaged exchange couplings J\ (shaded cir- 
cles) and J2 (open circles). Middle panel: convergence of the 
HTSE, measured as the ratio of the ninth-order term to the to- 
tal susceptibility. Bottom panel: saturation fields for the solu- 
tions with Ji < 0, J2 > (open triangles) and Ji > 0, J2 < 
(filled green triangles); the shaded stripe shows the experi- 
mental value /j,oH 3 ~ 23.4 T with the error bar of 0.2 T (see 
Ref. [24] for details) . The upper and middle panels refer only 
to the solution (a) with Ji < 0, J2 > 0, while the bottom 
panel shows both the solutions. 



data at sufficiently high temperatures, because the high- 
temperature limit constrains additional variable param- 
eters, such as xo and 9- I n the following, we will show 
that the specific heat data can also be a key for the as- 
signment of the Ji and Ji values, although such data are 
less sensitive to weak changes in the exchange couplings. 
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FIG. 5: (Color online) Specific heat of PbZnVO(P0 4 ) 2 : mag- 
netic contribution (open circles) and HTSE predictions for 
different solutions according to Table Mil The inset shows 
the experimental data (open triangles) along with the phonon 
(dashed line) and magnetic (solid line) contributions. 



B. Specific heat 

The temperature dependence of the specific heat (C p ) 
for PbZnVO(P04)2 is shown in the inset of Fig. [5] 
At low temperatures, C' p increases and has a peak at 
Tm = 3.9 K. This peak matches the bends of the suscep- 
tibility curves (Fig. [3]) and indicates a magnetic phase 
transition. Above Tjv, specific heat increases and shows 
a pronounced bend at 7.5 K. This region can be inter- 
preted as a sum of the rapidly increasing phonon con- 
tribution (Cphon) and the magnetic contribution (C mag ) 
that shows a maximum due to the correlated spin ex- 
citations right above Tjy. Since the exchange couplings 
in PbZnVO(PC>4)2 are relatively weak, one can expect 
C mag to reach its high-temperature limit (~ 1/T 2 ) above 
18 — 20 K. Then, the data can be fitted with the expres- 
sion 

(3) 

where the first term describes the high-temperature part 
of C mag , while further terms are Debye functions with 
Debye temperatures 9^ and integer coefficients c„. A 
similar model for C p hon has been used in our previous 
studies of vanadium phosphatesj 18 ' 19 i 41 ' 42 

Equation (JU) leads to a perfect fit of the experimental 
data above 18 K. The resulting A = 353(41) J K/mol 
has a large error bar, because at high temperatures C mag 
gives a small contribution to the total specific heat. The 
A value is in reasonable agreement with the J values: 
according to Ref. [llj, A/T 2 is the leading term in the 



HTSE for the specific heat, and 

A=(3R/8)(J 2 1+ Jl). (4) 

The solution (a') yields A ~ 396 J K/mol, while the 
solutions (b') and \b") lead to A = 220 - 240 J K/mol 
which is well below the experimental value. However, the 
distinction between the solutions (a 1 ) and (a") is hardly 
possible, because both the respective A values (396 and 
314 J K/mol) calculated from Eq. ((4|) agree reasonably 
well to the fitted value of 353 J K/mol. 

(n) 

The 9 D values are further used to estimate C p hon in 
the low-temperature region and to subtract the phonon 
contribution from the total specific heat. The resulting 
magnetic contribution is shown in Fig. [5l To test the 
correctness of the procedure, we estimate magnetic en- 
tropy S by integrating the C mag /T curve. We obtain 
S ~ 5.04 J/mol K, i.e., the magnetic entropy is slightly 
underestimated with respect to the expected value of 
i?ln2 ~ 5.76 J/mol K. This underestimate amounts to 
12 % and can be ascribed to a slight overestimate of 
the low-temperature phonon contribution and to a small 
amount of an impurity phase. 

The extracted experimental C mag (T) curve can be 
compared to the HTSE predictions^ 1 - At high temper- 
atures, the error bar for C mag is very high, because the 
measured heat capacity is dominated by the phonon con- 
tribution. Therefore, the direct fit of the C mag (T) curve 
is impossible. Nevertheless, one can use the intermediate 
region (10 — 20 K) and compare the data with the results 
of the susceptibility fitting. In Fig. El we plot the HTSE 
for the solutions (a') and (&') (Table IIIip . The curve 
for the solution (a') is in reasonable agreement with the 
experimental data (the agreement is further improved if 
the entropy weight error of about 10 % is corrected). The 
solution (&') clearly overestimates C mag below 15 K. Al- 
though both the solutions are obtained by fitting the sus- 
ceptibility data above 16 K, their convergence domains 
extend to lower temperatures (see the middle panel of 
Fig. 2] for the convergence of the susceptibility series). 
Therefore, the comparison in the region around 15 K is 
still reasonable. Due to the limited temperature range of 
the applicable specific heat data, distinguishing between 
(a') and (a") is hardly possible. 

Below Tjv, the specific heat of PbZnVO(PC>4)2 rapidly 
decreases. Such a behavior is qualitatively consistent 
with antifcrromagnetic ordering. However, the tempera- 
ture dependence of C mag in the ordered phase does not 
follow a simple T 3 trend expected for a regular antiferro- 
magnet. The origin of this behavior is presently unclear. 
Similar deviations from the T 3 behavior were observed 
in the other FSL compounds^ In these compounds, the 
columnar AFM ordering is stabilized by quantum fluctu- 
ations that select the collinear ground state among nu- 
merous ground states, which are degenerate in the classi- 
cal Heisenberg model (the so-called order-from-disorder 
mechanism) £ One might suggest that this feature mod- 
ifies the magnetic excitation spectrum and changes the 
low-temperature behavior of the specific heat. To clarify 
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this issue, further experimental and theoretical studies 
are desirable. 

To conclude the experimental section, the mag- 
netic susceptibility and the specific heat data for 
PbZnVO(PC>4)2 can be interpreted within the FSL 
model. Thermodynamic properties consistently point 
to ferromagnetic NN and antiferromagnetic NNN cou- 
plings in PbZnVO(PC>4)2 with the accurate estimates of 
Ji and J2. The reference to the saturation field sug- 
gests the spatial anisotropy of the NN couplings. In the 
following sections, we will use a microscopic approach 
in order to get further insight into the spin system of 
PbZnVO(PC>4)2 and to elucidate the relation between 
the structural features of this compound and individual 
exchange couplings. 

V. EVALUATION OF THE EXCHANGE 
INTEGRALS 

The calculated band structure of the material can be 
used to evaluate individual exchange couplings and to 
construct a microscopic model of magnetic interactions. 
In the case of magnetic insulators, there are two well- 
known approaches for the evaluation of the exchange 
couplings. The first approach utilizes the uncorrelated 
(LDA) band structure which is further mapped onto a 
TB model and, subsequently, onto a Hubbard model 
in the strongly correlated regime. Then, at half-filling 
the low-lying excitations of the system are properly de- 
scribed within a Heisenberg model, while the parame- 
ters of this model are expressed via microscopic param- 
eters of the starting Hubbard model. The alternative 
approach is based on a correlated band structure which 
is usually obtained via local spin density approxima- 
tion (LSDA)+£7 calculations that treat correlations in 
a mean-field way. To estimate exchange couplings, one 
can calculate total energies for ordered spin configura- 
tions and map the resulting energies onto the classical 
Heisenberg model. However, alternative methods, utiliz- 
ing Green's functions 4 ^ and a more sophisticated treat- 
ment of electronic correlations^ are also possible. 

The key feature of the second approach is the sizable 
and implicit dependence of the resulting exchange cou- 
plings on the parameters describing electronic correla- 
tions - in particular, on the on-site Coulomb repulsion 
parameter Ud- While in many systems the exchange 
couplings are sufficiently large and the uncertainty in 
Ud leads to 10 — 20 % uncertainty in J, the materials 
with weak exchange couplings are more difficult to study. 
In previous works ; 21 ' 42 we have shown that a change of 
1 eV in Ud can modify the signs of the exchange cou- 
pling constants in vanadium phosphates, thus different 
magnetic ground states are obtained. A careful fitting of 
the Ud value to the experimental data (e.g., the Curie- 
Weiss temperature) helps to establish the correct sce- 
nario. Still, it is preferrable to use the model approach 
and to keep the explicit dependence on the Coulomb re- 
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FIG. 6: (Color online) LDA density of states plot for 
PbZnVO(P0 4 )2: filled regions denote the total DOS, V and 
Zn contributions. The contributions of Pb, O, and P are 
shown by dark solid, light solid, and dashed lines, respec- 
tively. The Fermi level is at zero energy. 

pulsion parameter, while studying materials with weak 
exchange couplings. Another advantage of the model ap- 
proach is the insight into the microscopic mechanism of 
the exchange coupling, because contributions of different 
orbitals are easily separated. 

Following the model approach, we consider the uncor- 
related (LDA) band structure. The calculated density of 
states (DOS) is shown in Fig. [6] The valence bands of 
PbZnVO(P04)2 are formed by oxygen 2p orbitals with 
sizable contributions from Pb 6s states at —10 eV and 
from Zn 3d states at about —7 eV. The narrow bands 
at the Fermi level have vanadium id xy origin, while the 
other isolated band complex at 0.5 — 0.9 eV is formed 
by vanadium 3d xz and 3d yz states (see also Fig. [7|) . The 
e g states of vanadium lie above 1 eV and show strong 
hybridization with the 6p states of lead. The obtained 
crystal field splitting of the vanadium orbitals is consis- 
tent with the square-pyramidal (or distorted-octahedral) 
environment. Similar splittings have been observed in 
other vanadium compounds ! 21 ' 

The minimal microscopic model should include 8 d xy 
bands lying near the Fermi level (Fig. [7]). These bands 
originate from 8 vanadium atoms in the unit cell and 
are close to double degeneracy due to the weak inter- 
layer couplings. The d xy bands are fitted with a one- 
orbital TB model. The hoppings fv-^v Q f this TB 
model (Table IIV[) are further introduced into a one-band 
Hubbard model with the effective on-site Coulomb re- 
pulsion U c s- The [/off parameter is of the order of sev- 
eral eV . 21 ' 45 Then, the t <C U e s condition and the half- 
filling regime justify the mapping onto the Heisenberg 
model for low- lying excitations, and AFM contributions 
to the exchange integrals are estimated via the well- 
known expression J^ FM = At^j/U e g. 

To account for FM contributions to the exchange in- 



TABLE IV: Parameters of the TB model and the resulting 
exchange couplings in PbZnVO(P04)2. t xy ^ xy denote the 
hoppings between d xy orbitals, J AFM and J FM are antifer- 
romagnetic and ferromagnetic contributions to the total ex- 
change J, respectively. 







jrAFM (R) 


J FM (K) 


J(K) 


Jl 


11 


1.3 


-4.6 


-3.3 


J'x 


14 


2.0 


-4.2 


-2.2 


Jl 


-46 


21.9 


-0.1 


21.8 


Jk 


-38 


15.0 


-0.1 


14.9 



tegrals, one has to extend the model by the inclusion of 
unoccupied states. The underlying physical reason is as 
follows. If an electron hops to a half-filled orbital, its 
spin should be antiparallel to the spin of the electron on 
the destination site, hence AFM coupling is realized. In 
contrast, the hopping to an empty orbital favors parallel 
arrangement of spins due to the Hund's coupling (on-site 
exchange) on the destination site. Analytical expressions 
for the FM coupling of this type were first derived by 
Kugel and Khomskii as a special case of their model. 46 
Later on, this expression was employed by Mazurenko 
et a/4£ for the analysis of the exchange couplings in 
Na 2 V 3 07. Following Refs. |H and H3, we estimate the 
FM contributions to the exchange couplings as 

t-FM = _V^ 1 jj > 

ij ^ (C/eff + A Q ) (U cS — J c ff + A Q ) ' [ ' 

where the index a indicates empty orbitals, t x ^ a de- 
notes the hopping from the d xy orbital on site i to an 
a orbital on site A a = e a — e xy is the energy split- 
ting between the occupied (xy) and empty (a) orbitals, 
U e g and J e g are the effective on-site Coulomb repulsion 
and the on-site (Hund's) exchange parameters, respec- 
tively. Hereby, we assume that the on-site repulsion and 
exchange are similar for the xy and for all the a orbitals. 

The hoppings t xy ^ a are estimated from a fit of the 
LDA band structure with a multi-orbital TB model. For 
the evaluation of J FM , we include in the model both the 
d yz and d xz states of vanadium (see Fig. [7|) and omit 
the e g states due to their strong hybridization with lead 
orbitals. This approach is justified by a reference calcu- 
lation for the isostructural SrZnVO(P04)2 (not shown), 
where the e g states are weakly hybridized and show small 
hoppings to the xy states (t xy ~* e <> < 5 meV). Another 
justification is given by an extended TB model includ- 
ing all the vanadium 3d and, additionally, lead 6p states. 
Due to the presence of uncorrelated (Pb 6p) states, this 
extended model can not be directly mapped onto the 
Heisenberg model, and one has to reduce the model to 
five vanadium orbitals first. Still, the parameters of the 
extended model suggest weak hoppings between the half- 
filled d xy orbitals and the e g orbitals of vanadium or the 
6p orbitals of lead (t xy ^ e <> , 1fv^ ph6 P < 2 meV). There- 
fore, the xy — y e g hoppings in the reduced (five-orbital) 
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FIG. 7: (Color online) LDA band structure of 
PbZnVO(P0 4 ) 2 (thin light lines) and the fit of the 
three-orbital TB model (thick lines). The Fermi level is at 
zero energy. The notation of k points is as follows: r(0, 0, 0), 
X(0.5,0,0), M(0.5, 0.5,0), Y(0, 0.5,0), Z(0, 0,0.5), and 
T(0.5,0,0), where the coordinates are given in units of the 
respective reciprocal lattice parameters. Note that most of 
the bands are close to double degeneracy due to the weak 
interlayer coupling. 



model should also be weak. 

For the U e s and J c s parameters, we use the represen- 
tative values of 4.5 eV and 1 eV, as employed in the pre- 
vious model analysisi 21 i 42 ' 45 i 47 We should note that the 
J e ff value has little influence on J FM , and the main un- 
certainty arises from U e s (J FM roughly scales as 1/U^ S ). 
However, the change of U e g in the reasonable range of 
4 — 6 eV leads to a moderate change in the resulting J 
values, despite all the couplings are weak. Such a moder- 
ate dependence on U s is a clear advantage of the model 
approach as compared to the LSDA+C/ calculations (see 
Ref.HH). 

The calculated exchange integrals arc listed in Table HVl 
along with their FM and AFM contributions. We find 
that the NN couplings J\ and J[ are FM, while the 
NNN couplings J2 and J' 2 are AFM in remarkable agree- 
ment with the experimental results: J\ ~ —5.2 K and 
J 2 — 10.0 K. Moreover, even the difference J{ — J\ ~ 
1.1 K perfectly matches the experimental value of 1.1 K. 
The interactions J\ and J{ show both FM and AFM 
contributions, although the former contribution domi- 
nates. In the case of J2 and J 2 , the values are some- 
what overestimated as compared to J2 ~ 10 K. This 
effect is also observed for the other AA'VO(P04)2 vana- 
dium phosphates^ The NNN interactions show negligi- 
ble FM contributions consistent with long V-V separa- 
tions of about 6.5 A. The largest interlayer coupling is of 
the order of 0.01 K. 

Apart from the support of the experimental results, 
our microscopic model provides further insight into the 
spin system of PbZnVO(P04)2- In particular, the com- 
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putational study gives independent estimates of J 2 and 
J' 2 . We find pronounced spatial anisotropy of the NNN 
couplings (J2/J2 — 0.68) which is comparable to that 
in BaZnVO(P0 4 ) 2 (J2/J2 ^ 0.61) and Pb 2 VO(P0 4 ) 2 
(J2/J2 ^ 0.67)^- In Ref. EH, we have shown that the 
spatial anisotropy of the NNN couplings is mainly deter- 
mined by the difference in V-0 distances for the respec- 
tive superexchange pathways. In the next section, we will 
discuss in more detail the geometry of the superexchange 
pathways in PbZnVO(P04)2 and the relation between 
the structural features and magnetic interactions in FSL 
compounds. 



VI. MODELING 

The magnitudes of superexchange interactions depend 
delicately on geometrical parameters of the superex- 
change pathway. In the conventional M-O-M superex- 
change scenario (M is a transition metal), the exchange 
integral is most sensitive to the angle at the oxygen 
atom. 49 However, more complicated superexchange path- 
ways can show different trends. In the case of vanadium 
phosphates (and, more generally, any transition-metal 
compounds with XO4 tetrahedra, X being a main-group 
cation), the typical scenario is the "superexchange via a 
single tetrahedron" , as shown in the left panel of Fig. [5] 
This scenario corresponds to the interactions J 2 and J' 2 
in the layered vanadium phosphates AA'VO(PC>4) 2 . Ac- 
cording to Ref. Hi], the A and A' metal cations have influ- 
ence on the: (i) buckling of the [VOPO4] layer, i.e., on the 
angles at oxygen atoms that determine the orientation of 
the tetrahedron with respect to the VO5 pyramid; (ii) 
stretching of the [VOPO4] layer, i.e., vanadium-oxygen 
distances. Our microscopic analysis has shown that the 
exchange integrals are more sensitive to the latter, while 
the change in the angles (in a reasonable range) leads to 
a weak modification of the exchange couplings (see Fig. 8 
in Ref. 2li). These findings suggest that the 0-(P)-0 in- 
teraction is rather insensitive to the mutual orientation 
of the VO5 pyramids and the PO4 tetrahedron. How- 
ever, the microscopic origin of this interaction remains 
unclear. 

We start with the careful analysis of individual struc- 
tural parameters of PbZnVO(PC>4) 2 . To characterize the 
J 2 and J'2 superexchange pathways, we average the re- 
spective V-0 distances (Table |rfj}: d ~ 2.01 A (J 2 ), 
d! ~ 2.05 A (4). According to Ref. EH, these num- 
bers mainly determine the magnitudes of J 2 and J^. 
The d value is somewhat higher as compared to the 
other AA'VO(P0 4 ) 2 compounds [e.g., d ~ 1.98 A in 
SrZnVO(P04) 2 ]. In consequence, the interaction J 2 
in PbZnVO(P04) 2 should be weaker than in the SrZn 
compound. However, this is not the case: we find 
J 2 ~ 22 K (Table HVj) in contrast to J^ FM ~ 19 K in 
SrZnVO(P04) 2 . 21 Thus, other structural factors influ- 
ence on the J 2 value. Indeed, several P-0 distances in 
PbZnVO(P04) 2 are smaller than in the other compounds 



of the family. This gives the clue that the magnitude of 
the superexchange interaction is also controlled by the 
structural parameters of the PO4 tetrahedron. 

To analyze the issue in more detail, we perform calcu- 
lations for model structures. Such structures include the 
magnetic [VOPO4] layers separated by Li cations. The 
details of the procedure can be found in Ref. l2l|. We 
assume that the magnetic interactions are determined 
by the structure of the [VOPO4] layer, while metallic 
cations simply provide proper charges (then, the com- 
plex [AA'P04] +1 block can be substituted by a layer of 
Li +1 cations). This assumption perfectly holds for the 
AFM interactions (i.e., for the xy — > xy hoppings), 21 
while the case of the FM interactions is more problem- 
atic (hoppings between different d orbitals depend on the 
structure of interlayer blocks). In the following, we fo- 
cus on the AFM interactions (J 2 and J2) only. Clearly, 
the FM contributions to these interactions are negligible 
due to the large V-V separation (see Table IIV|) , hence 

First, we consider the uniform scaling of individual 
polyhedra in the magnetic layer and trace the change in 
the coupling constant J 2 (due to the uniform scaling, the 
layer remains regular, hence J 2 = J^)- Both the pyramid 
and the tetrahedron exhibit sizable influence on the ex- 
change integral (see middle panel of Fig. [8]). The trend for 
the pyramid basically reproduces the change in J'^ in the 
AA'VO(P04) 2 compounds: d! is increased from 2.00 A 
in BaCdVO(P0 4 ) 2 up to 2.08 A in SrZnVO(P0 4 ) 2 ; J' 2 is 
correspondingly reduced from 17 K down to 7 K.- 2 * The 
curve for the tetrahedron shows a smaller slope. Thus, 
the exchange coupling is less sensitive to the intermediate 
0-(P)-0 link, and the V-0 distances (i.e., the overlap 
between vanadium and oxygen orbitals) are more impor- 
tant. 

Now, we will elucidate the relevant geometrical param- 
eter of the tetrahedral group. Uniform scaling of the 
tetrahedron implies that both the P-0 distances and the 
0-0 distances (edges of the tetrahedron) are increased. 
In the next run of model calculations, we kept the 0-0 
distances fixed and varied the P-0 distances by shift- 
ing the phosphorous atom inside the tetrahedron. The 
shift of the P atom leads to the non-equivalence of J 2 
and J 2- However, the J2/J2 remains close to unity even 
if the difference between the respective P-0 distances is 
pronounced (see the right panel of Fig. [5]). In contrast, 
the change in individual V-0 distances leads to the rapid 
decrease of J' 2 / ' J%. This result enables us to establish the 
V-0 distances in the pyramid and the 0-0 distances in 
the tetrahedron as the relevant geometrical parameters 
for the superexchange interactions J 2 and J^- 

A more general result of the above analysis deals with 
the contribution of the phosphorous atom to the superex- 
change. Since the P-0 distances have little influence on 
the exchange integrals, we can conclude that the phos- 
phorous orbitals play a minor role in the coupling. The 
superexchange is the result of the interaction between the 
oxygen atoms on the edge of the tetrahedron. This pic- 
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FIG. 8: (Color online) The superexchange pathway comprising the single tetrahedral PO4 group (left panel) and the evolution 
of the exchange couplings upon variation of the geometry of this pathway (middle and right panels). The middle panel shows 
the change in the exchange coupling J2 along the scaling of the polyhedra in a uniform way, while the right panel presents the 
spatial anisotropy of the exchange couplings (J2/J2) due to the different P-0 or V-0 distances in the respective polyhedra 
(Ad is the difference between the distances, relevant for J2 and J^)-— While the parameters of one polyhedron (pyramid or 
tetrahedron) are varied, the geometry of the other polyhedron is fixed at d(V-O) = 1.95 A for the pyramid and d(P-O) = 1.54 A 
for the tetrahedron (the crossing point of the curves in the middle panel). 
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FIG. 9: (Color online) Plot of the Wannier functions for the 
superexchange pathway, including two VO5 pyramids and the 
single PO4 tetrahedron. Arrows indicate the overlap of oxy- 
gen p orbitals. 



ture is illustrated by the plot of the Wannier functions in 
Fig. [5J Each Wannier function is composed of the vana- 
dium d xy orbital and oxygen p orbitals. As expected, 
the contribution of phosphorous orbitals is minor. Then, 
the superexchange can be considered as the 7r-overlap of 
vanadium and oxygen orbitals with a further 7r-overlap of 
the two oxygen orbitals on neighboring atoms. This sce- 
nario also accounts for different changes in J 2 upon scal- 
ing the pyramids and the tetrahedra (see middle panel of 
Fig. [5]). The scaling of the pyramids modifies two V-0 
distances and has stronger influence on the exchange as 
compared to the scaling of the tetrahedra that leads to 
the change in a single 0-0 distance. In the next section, 
we will discuss the implications of this result for a tuning 
of the magnetic interactions in FSL materials. 



VII. DISCUSSION AND SUMMARY 

Our studies identified PbZnVO(P04)2 as a FSL ma- 
terial with ferromagnetic NN interactions, antiferromag- 
netic NNN interactions, and a spatial anisotropy of the 
exchange couplings. Thermodynamic measurements and 
band structure calculations provide convincing evidence 
for this scenario and establish solid and accurate esti- 
mates of individual exchange couplings. Now, one can 
compare PbZnVO(P04)2 to other FSL materials and ob- 
tain deeper experimental insight into the properties of the 
spatially anisotropic FSL model. 

Table [V] presents key characteristics of the magnetic 
susceptibility and the specific heat for a number of FSL 
materials with FM J\ and AFM J 2 ■ These materials re- 
veal different magnitude of the frustration, as evidenced 
by the effective frustration ratio a = J%j J\. The a value 
is reduced from BaCdVO(P0 4 ) 2 to SrZnVO(P0 4 ) 2 and, 
further, to Pb 2 VO(P0 4 ) 2 and PbZnVO(P0 4 ) 2 . The 
two latter compounds show similar a ~ —1.9. How- 
ever, they are different with respect to their NN cou- 
plings. The saturation field measurement indicates 
that Pb 2 VO(P0 4 ) 2 is spatially isotropic (J 4 ~ J[), while 
in the case of PbZnVO(P0 4 ) 2 J[ - Ji ~ 1 K. Neverthe- 
less, the positions and magnitudes of the susceptibility 
maxima in the two compounds are close to each other. 
On the other hand, the enhancement of the frustration in 
SrZnVO(P0 4 ) 2 and BaCdVO(P0 4 ) 2 leads to the shift of 
the maximum to lower temperatures and to the increase 
of the susceptibility value at the maximum. 

At first glance, the characteristics of the specific heat 
are slightly different: the positions of the maxima show 
the same trend as observed in the susceptibility data. 
Yet the maximum specific heat value in PbZnVO(P0 4 ) 2 
is lower than that in Pb 2 VO(P0 4 ) 2 and comparable to 
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TABLE V: Basic characteristics of thermodynamic data for the FSL compounds with FM Ji and AFM J 2 . a = J2/J1 is the 
effective frustration ratio, T^ax and T^ &K are the positions of the magnetic susceptibility and the specific heat maxima, while 
X max and (7 max are the absolute values at the maxima. J c = + Ji) 1 ^ 2 defines the thermodynamic energy scale. 



Compound 


a 


rpX IT 
1 raax/ ° c 


X max J c 


T° 1 1 

± max/ °c 


C max /i? 


Ref. 


PbZnVO(P0 4 ) 2 


-1.9 


0.78 


0.172" 


0.44 


0.41" 




Pb 2 VO(P0 4 ) 2 


-1.9 


0.83 


0.189 


0.42 


0.45 


17 and .18 


SrZnVO(P0 4 ) 2 


-1.1 


0.52 


0.303 


0.28 


0.41 


18 


BaCdVO(P0 4 ) 2 


-0.9 


0.52 


0.346 


0.31 


0.31 


19 



"These values are reduced for about 10 % due to diamagnctic 
impurities, see text for details. 



C max in SrZnVO(P0 4 ) 2 . Still, this result is likely in- 
significant, because the C m&x value is affected by the im- 
purities. In Sec. IIV1 we have shown that both the mag- 
netization and the specific heat data point to ~ 10 % of 
diamagnetic impurities in the samples under investiga- 
tion. Once this error is corrected, we find C max ~ 0.45i? 
which coincides with the result for Pb2VO(P0 4 )2- More- 
over, x max ^c rescales to 0.191 in perfect agreement with 
x max J c = 0.189 in Pb 2 VO(P0 4 ) 2 . 

Our comparison of the different FSL materials indi- 
cates that the spatial anisotropy of the spin lattice has a 
minor effect on thermodynamic properties. The behavior 
of the magnetization and the specific heat is mainly de- 
termined by the frustration ratio that controls the mag- 
nitude of quantum fluctuations and determines spin cor- 
relations in the system. This finding provides experimen- 
tal verification for simulation results which consistently 
showed weak changes in thermodynamic properties upon 
the distortion of the FSL^ Thus, the distortion of the 
spin lattice is a secondary effect compared to the frus- 
tration. In the case of the FSL, low-symmetry materials 
with complex structures can be reliably considered as ex- 
perimental realizations of the regular model. 

Extending the above statements to the ground state 
properties and using a ~ —1.9, one would suggest that 
PbZnVO(P0 4 )2 undergoes long-range ordering towards 
columnar antifcrromagnetic state at TV — 3.9 K (in zero 
field). This suggestion is further confirmed by the ref- 
erence to Pb2VO(P0 4 )2 with similar exchange couplings 
(Ji ~ -5.1 K, J 2 ~ 9.4 K), similar T N ~ 3.5 K, and sim- 
ilar transition anomalies in the magnetic susceptibility 
and the specific heat i 17 ' 18 The columnar AFM ordering in 
Pb2VO(P0 4 )2 was directly confirmed by neutron scatter- 
ing and nuclear magnetic resonance measurements 1 51 1 52 
Thus, a similar ordering in PbZnVO(P0 4 )2 looks likely. 
Then, the role of the spatial anisotropy is mainly re- 
stricted to determine the direction of columns. Since 
\J[ \ < \Ji\, we expect that spins will show parallel align- 
ment along the b direction (interaction J\). To test this 
prediction, further neutron scattering studies are desir- 
able. We also hope that experimental realizations of the 
spatially anisotropic FSL will stimulate theoretical inves- 
tigation of the respective model. 

Finally, we turn to the structural aspects of the 



study. Different geometries of the [VOP0 4 ] layers in 
the AA'VO(P0 4 )2 compounds stimulated us to study 
the geometrical parameters that influence the exchange 
couplings. For the NNN interactions J2 and J' 2l these 
parameters are the V-0 and 0-0 distances. Mutual ori- 
entation of the VO5 pyramids and the P0 4 tetrahedron 
has minor effect on the exchange integrals (within a rea- 
sonable range of possible geometries). The role of the 
tetrahedrally coordinated X cation is to define the size of 
the tetrahedron and the relevant 0-0 distance. As the 
X cation gets larger, the NNN interactions in the FSL 
materials are reduced. This result helps to explain the 
change of J2 in the family of the FSL compounds. 

Vanadium phosphates normally show J2 of 9 — 
10 K 17 ' 18 i 19 in agreement with the small size of the 
P +5 cation (ionic radius r = 0.17 A, 34 typical O- 
O separations d Q -o = 2.3-2.5 A)£L2L22. As we 
turn to Li2VOSi0 4 with the larger Si +4 cation (r = 
0.26 A^i do-o = 2.63 AS), J 2 is reduced to ~ 
6 K. 1 " -18 In Li 2 VOGc0 4 , one finds an even smaller J 2 ~ 
4 KlSii 8 - due to the larger Ge+ 4 cation (r = 0.39 A^i 
d _ = 2.78 A 53 ). We should note that VOMo0 4 
with J2 ~ 20 K 13 ' 14 does not follow this trend due to 
the slightly different crystal structure and the differ- 
ent superexchange scenario in this compound. Mo +6 
is a transition-metal cation, hence its 4d orbitals have 
low energies and show stronger overlap with the vana- 
dium orbitals as compared to the p orbitals of the 
main-group cations (P, Si, Ge)^ Another exception is 
BaCdVO(P0 4 ) 2 with its unusually low J 2 ~ 3.2 K. The 
origin of this anomaly remains unclear and should be a 
subject of future investigations. 

Based on the above analysis, we can propose a new 
route towards strongly frustrated square lattice materi- 
als. In the AA'VO(P0 4 ) 2 , Li 2 VOSi0 4 , and Li 2 VOGe0 4 
compounds, the |Ji| < J2 regime is realized. To reach 
the critical regions at J2/IJ1I — 0.5, one has to reduce 
J2 (irrespective of the sign of J\). Since the J2 value is 
controlled by the size of the tetrahedron, one should use 
a larger tetrahedrally coordinated cation. The promis- 
ing candidate is As +5 with r = 0.34 A 3 -! and typical 
do-o — 2.65 — 2.80 A, providing J 2 of about 5 K or 
even lower. Our preliminary band structure calcula- 
tions for vanadium arsenates show weak contribution of 
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As states at the Fermi level, hence the magnitude of 
superexchange interactions will mainly depend on the 
0-0 distances. The chemistry of vanadium(IV) arsen- 
ates remains weakly explored (in particular, none of the 
AA'VO(As04)2 compounds have been reported) and de- 
mands further investigation. 

In summary, we have reported the crystal structure to- 
gether with an experimental and microscopic study of the 
magnetic properties of PbZnV0(P04)2- This compound 
can be understood as a frustrated square lattice mate- 
rial with Ji ~ —5.2 K, J 2 ~ 10.0 K, and sizable spatial 
anisotropy of nearest-neighbor couplings (J[/Ji — 0.8) 
and next-nearest-neighbor couplings (J2/J2 — 0.7). The 
spatial anisotropy shows little effect on thermodynamic 
properties. At T N ~ 3.9 K, PbZnVO(P0 4 ) 2 under- 
goes a phase transition which is likely associated with 
columnar antiferromagnetic ordering, similar to the pre- 
dictions of the regular model for the frustration ratio 
a = J2/J1 — —1.9. Our microscopic study helps to 



identify the relevant geometrical parameters for the next- 
nearest-neighbor interactions in FSL materials. We show 
that the magnitudes of such interactions crucially depend 
on the size of the tetrahedrally coordinated cation. This 
finding suggests a new route towards strongly frustrated 
square lattice materials in the critical region of the re- 
spective phase diagram. 
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